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Abstract 

We estimate the entropy associated to a background of squeezed cosmic gravitons, 
and we argue that the process of cosmo logical pair production from the vacuum 
may explain the large amount of entropy of our present universe. 
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Entropy production in the cosmological amplification of the vacuum 
fluctuations 

A successful inflationary model is expected to solve the various problems of 
the standard cosmological scenario and then, in particular, must provide some 
explanation to the unnaturally large value of the entropy to-day stored in the 
cosmic microwave background (CMB). The solution to the entropy problem relies, 
conventionally, on the so-called "reheating" process [1], occurring at the end of 
the inflationary era, and realized in various model-dependent ways such as bubble 
collisions and/or inflaton decay. In view of the difficulties encountered by some 
model in providing an efficient enough reheating phase [1,2], however, it seems 
interesting to point out an additional source of entropy in the inflationary scenario: 
the parametric amplification of the vacuum fluctuations, or (equivalently stated, 
in a second quantization formalism) the pair production from the vacuum, due to 
changes in the background metric field. 

As discussed in previous papers [3-5], the growth of the average number of 
field quanta, due to pair production, is naturally associated to the growth of 
the entropy of that field. On the other hand the pairs spontaneously emitted 
from the vacuum, in an evolving cosmological background, necessarily appear in a 
"squeezed" quantum state [6,7], namely in a state in which the quantum fluctua- 
tions of a given field operator Q are suppressed with respect to the vacuum, while 
the fluctuations of the canonically conjugate operator P (called " superfluctuant" ) 
are correspondingly enhanced (see, for instance, [8,9]). 

In order to compute the entropy associated to cosmological pair production, 
we shall thus propose in this paper a " coarse graining" approach to non-equilibrium 
entropy, valid for squeezed states, in which the loss of information associated to 
the reduced density matrix is represented by the increased dispersion of the su- 
perfluctuant operator P. Such a quantum-mechanical approach is different from 
a recent field-theoretic computation of the entropy of the cosmological pertur- 
bations [10,11]. The two methods provide results which coincide in the large 
squeezing limit; our expression for the entropy, however, is always positive definite 
(unlike the classical approach of [10]), and seems then to be valid even in the small 
squeezing regime. 

The definition of entropy we shall provide holds in general for any mechanism 
of particle production from the vacuum (more generally, from any coherent ini- 
tial state), not necessarily of gravitational origin, which can be represented as a 
mixing of positive and negative frequency modes through an appropriate Bogoli- 
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ubov transformation. We shall apply this definition, in particular, to estimate the 
entropy change associated to the formation of a background of cosmic gravitons. 
By using a very general model of cosmological evolution, we shall find an entropy 
which may be as large as the CMB electromagnetic entropy, provided the transi- 
tion from inflation to the radiation-dominated phase occurs at a curvature scale 
of the Planck order. 

Let us recall, first of all, that the quantum description [6,7] of the amplifi- 
cation of scalar or tensor fluctuations is based on the separation of the field into 
background solution and first order perturbations, and on the expansion of the so- 
lution to the perturbed wave equation into \in > and \out > modes. The complex 
coefficients of this expansion are interpreted, in a second-quantization formalism, 
as annihilation and creation operators for a particle (6, &t) and the corresponding 
anti-particle (b,W). The relation between \in > and \out > mode solutions can 
thus be expressed, for each mode k, as a Bogoliubov transformation between the 
\in > operators (6, b\ b, bt) and the \out > ones (a, a\ a, a)): 

a k = c+(k)b k + c*_(k)bi k 
al k =c-(k)b k +c* + (k)~bl k 

where |c + | 2 — |c_| 2 = 1. By parametrizing the Bogoliubov coefficients c± in terms 
of the two real numbers r > and 6>, 

c+(k) = coshr(k) , c*_ (k) = e 2l6k sinhr(/c) (2) 

the relations (1) can be re- written as unitary transformations generated by the 
(momentum-conserving) two- mode squeezing operator 

E fc = eM4hb- k - z k b{bl k ) , z k = r{k)e 2i ° k (3) 

(r is the so-called squeezing parameter) as 

a k = T, k b k T,l (4) 

(and related expressions for a^,a^ and a). Starting with the \in > vacuum state 
|0 >, such that b k \0 >= = b- k \0 >, the pair production process leads then, for 
each mode, to the "squeezed vacuum" state [6,7] \z k >= E fc |0 >, which satisfies 
a k\ z k >= = d- k \z k >, and which has an averaged particle number 

n k =< z k \b\b k \z k >=< z k \tf_ k b- k \z k >=< 0|a[a fc |0 >= \c-(k)\ 2 = sinh 2 r(/c) 

(5) 
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For the purpose of this paper it is convenient to note that the two-mode trans- 
formation (4) can be factorized as the product of two one-mode transformations, 
as usually done in squeezed-state theory [9] , by defining a new pair of annihilation 
operators (ai,a\) through the complex rotation 

a = —j=(ai — wi) , a = —=(ai + iai) (6) 
V 2 V 2 

and a similar expression for 6, 6, 

b=^(b 1 -ib 1 ) , 6=-^(&! + i6i) (7) 

(obviously, by definition, [ai,a{] = [ai,a{] = [&i,&i] = [&i,&i] = 1, = 
[&i,&i] = 0). The transformation (4) can thus be rewritten, for each mode, as 

a = h coshr + b\e 2i9 sinhr - coshr + b\e 2te sinhr) = Ei&iE{ - zE^E} (8) 

where Ei and Ei are the one-mode squeezing operators for b± and b\, 

Si = exp(^6? - Z -b\ 2 ) , E 1 =exp(|6?-|6l 2 ) (9) 

Given any two-mode squeezed state \z >, with squeezing parameter r > 0, 
it is always possible to introduce two operators x and x whose variances (Ax) z , 
(Ax) z are amplified with respect to the vacuum value (Ax)o, (Ax)o, namely [8,9] 
(Ax) z /(Ax) = (Ax) z /(Ax) = exp(r), where (Ax 2 ) z =< z\x 2 \z >-(< z\x\z >) 2 
(and the same for x). In terms of these operators, b\ and b\ have the differential 
representation 

i{B+ir/2) ^ i(e+n/2) 

h= (x + d x ) , h= (x + ds) (10) 

where the phase has been chosen according to eq.(2), in such a way to identify the 
x and x operators with the superfluctuant ones. 

Therefore, in the (x, x)-space representation the wave function ifj Zk for the two- 
mode squeezed vacuum, determined by the condition a^ 2 = 0, can be factorized 
according to eq.(8) as 

?/> 2fe (x, x) =< xx\z k >= (x)ipl k (x) (11) 
where ip\(x) is fixed by the differential equation 

cosh r(x + dx)^ 1 = sinh r(x — dx)^ 1 (12) 
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With the normalization condition < Zk\zk >= 1 one obtains 

^ Zk ( x ,x) = (^exp[-^(x 2 +x*)] (13) 

7T Z 

where the real number a k = exp[— 2r(k)] measures the departure of ip Zk from the 
vacuum wave function (a = 1). 

The dynamical evolution leading to the formation of squeezed states is thus 
accompanied, in general, by a loss of information corresponding to a larger disper- 
sion of x and x around their mean values, Ax = Ax = 1/ \pla. Such an increase 
in uncertainty is measured, in the (x,x) representation, by the "flattening" of 
the gaussian (13) with respect to the vacuum wave function, with a probability 
distribution P k (x,x) = \ < xx\z k > | 2 which defines the reduced density operator 



Pk 



/,«;•,> " " : * " ; "''- r " " ' Trpk = X (14) 



The entropy production associated to the transition |0 >^ \z k >, according to 
this reduction scheme, can now be computed by applying the usual information- 
theoretic definition of entropy (S = — Trp In p), and by subtracting the constant 
contribution of the vacuum state, 

S = / dxdxe-( x2+i2) ln(- ) = 1 + hi7r (15) 

TV J 7T 

The result can be simply expressed in terms of the averaged particle number 

as 

S(k) = —Trpk In pk — So = — J dxdxPk(x, x) lnPk(x, x) — So = 

= 2r(k) = 2sinh" 1 |c_(fc)| = 2 ln(VW^ + y/1 + n k ) (16) 

This expression for the entropy density per mode holds even if the initial 
vacuum is replaced in general by a coherent state \a k >, such that bk\ak >= 
<*k\ a k >; moreover, it is always positive definite, and can be applied for any value 
of the squeezing parameter r. In the large (or, equivalently, large squeezing) 
limit, eq.(16) gives 

S(k) — In |c_(fc)| 2 = lnn fc (17) 

which coincides exactly ( for n k » 1) with the entropy associated to an oscillator 
in a state of thermal equilibrium, and with the result of [10,11] (such a relation 
between the entropy and the mixing Bogoliubov coefficient was first suggested 
in [3]). The same result can be obtained if the reduction procedure, from the 
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pure squeezed state down to a statistical mixture, is performed (mode by mode) 
in the Fock space spanned by the eigenstates \n k h- k > of the number operators 
N k = b\b k , N k = bl k b- k . 

Indeed, as stressed in [6], the number N is a superfluctuant operator in the 
state \z k >, while the variance of the conjugate phase operator is squeezed. For 
large values of r, the loss of information with respect to N is measured, as discussed 
also in [5,11], by the density operator 

Pk = ^2 Pk(n)\n k h- k >< n k h- k \ (18) 

n 

where P k (n) = \ < nh\z k > | 2 . It is well known, on the other hand, that 
any process of pair production from the vacuum described by the Bogoliubov 
transformation (1) is characterized by a probability distribution [12] 

\r (k)\ 2n 

*(») - ( i + Mil I')-*- < 19 > 

One thus obtain, for |c_| >> 1, that S(k) = —Trp k \rvp k ~ ln|c_(/c)| 2 , in agree- 
ment with eq.(17). 

The expression (16) provides the entropy density per mode for all cases of 
particle production (not necessarily of gravitational origin) parametrized by the 
Bogoliubov transformation (1). The total entropy S in a proper spatial volume V 
is thus obtained by summing the squeezing contributions over all modes of proper 
frequency uj (following the coarse graining approach of [4,5], we have neglected 
any correlation among modes) 



S=^ f 1 r{u)u 2 du (20) 

where ujq and u\ delimit the frequency interval in which the squeezing mechanism 
is effective. 

We shall apply this result in order to estimate the entropy growth associ- 
ated to the cosmological productions of gravitons in a homogeneous and isotropic 
background, whose scale factor a(t) describes the evolution from a primordial in- 
flationary phase to the standard radiation and matter-dominated expansion. In 
such case proper and comoving frequency (k) are related by uj = k/a(t), and the 
graviton spectrum is determined by the Bogoliubov coefficient as follows [13-15] 
(we follow the notations of [13]) 

\c-(u)\ ~ ( — )" l<s| , UJ 2 < UJ < UJi 

UJ i x , UJ n (21) 

|c_M|~( — )-l 5 l( — r 1 , UJ <UJ<UJ 2 
UJ\ U)2 
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Here S is an order of unity parameter depending on the kinematical behaviour 
of the inflationary phase (\5\ = 2, for instance, in the case of de Sitter infla- 
tion); u>i is the maximum cutoff frequency depending on the curvature scale Hi at 
the beginning {t = t\) of the radiation-dominated evolution; u>2 is the frequency 
corresponding to the radiation-matter transition {t = £2); finally, ~ Hq ~ 
10 _18 sec _1 is the minimum frequency fixed by the to-day value of the Hubble 
radius (uji >> UJ2 » uJq). 

In this example |c_| >> 1, so that r(oj) ~ ln|c_(a;)|. By integrating the 
squeezing parameter according to eq.(20) we find that the total entropy per co- 
moving volume is constant, just like the black- body entropy. By keeping only the 
dominant terms, we obtain for the graviton entropy S g in a comoving volume £ 3 

S 9 = (aM )3 (22) 

Neglecting numerical factors of order unity, the total entropy (aiuJi) 3 can be easily 
estimated at any given observation time to during the matter-dominated era, by 
rescaling down ui(to) to ujo as [13-15] ui(to)/uJo = ki/ko ~ (Hiai)/(Hoao), where 
Hi = H(ti), ai = a(ti) and so on. By recalling that a ~ if -1 / 2 in the radiation 
phase, we obtain 

Ul _ /#l\l/2/#202x _ Hi 1/2 MpOp 

~ { H 2 Kad{ H a )mat ~ { M P } { H a Q } { 6) 

where the transition curvature scale Hi has been conveniently expressed in units of 
Planck mass Mp. By introducing the black-body temperature T 7 , which evolves 
adiabatically (a(t)T 7 (t) = const) from T 7 (t ) to T 7 (tp) = Mp, we find from 
eq.(23) 

^i(to) = (^) 1/2 T 7 (to). (24) 

S g can thus be rewritten 

S 9 - H^-^) 372 (25) 

where «S 7 ~ [a(t)£T 7 (t)] 3 = const is the usual black-body entropy of the CMB 
radiation (in terms of the to-day parameters, (oo^T 70 ) 3 ~ (T l0 /H ) 3 ~ (10 29 ) 3 ). 

Eq.(25), which represent the main result of this paper, holds in general for 
particle production with a spectrum n(iv) ~ a; - ' 6 ', quite irrespective of the kind of 
fluctuations which are amplified (even for photons, in a suitable non-conformally 
flat background), and is only weakly dependent on the background kinematics 
(through the parameter \S\, which is typically of order unity [13-15]). Particle 
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production from the vacuum is thus a process able to explain the observed cos- 
mological level of entropy, provided the curvature scale at the inflation-radiation 
transition is of the order of the Planck one, Hi ~ Mp. 

In the standard inflationary scenario, based on an accelerated phase of de 
Sitter-like exponential expansion, there are well known phenomenological bounds 
[16] which imply Hi < 10~ 4 Mp, and which rule out the mechanism discussed here 
as a possible explanation of the observed entropy. However, as recently stressed 
in [13,17,18], such constraints are evaded if the de Sitter and/or the radiation- 
dominated phase are preceeded by a phase of accelerated evolution and growing 
curvature, like in the recently proposed " pre-big-bang" models [18,19] of a duality- 
symmetric string cosmology. We thus conclude that the result obtained in this pa- 
per provides an additional reason of interest for such models, where the maximum 
curvature scale is only constrained by [13,17,18] Hi < 10 2 M P . 
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Note added 

After this paper was submitted we received a preprint by R. Brandenberger, 
T.Prokopec and V.Mukhanov, "The entropy of the gravitational field" (Brown- 
HET-849, August 1992), whose content overlap to some extent with ours, and 
in which a result similar to our equation (25) is obtained through a different 
procedure. 
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